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Abstract 

We analyse numerically the onset of pre-big bang inflation in an inhomoge- 
neous, spherically symmetric Universe. Adding a small dilatonic perturbation 
to a trivial (Milne) background, we find that suitable regions of space undergo 
dilaton-driven inflation and quickly become spatially flat (fi — > 1). Numerical 
calculations are pushed close enough to the big bang singularity to allow cross 
checks against previously proposed analytic asymptotic solutions. 
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I. INTRODUCTION 



Increasing attention has recently been devoted to a possible alternative to the standard 
inflationary paradigm, the so-called pre-big bang (PBB) scenario [p]], 0, ||. While referring 
the interested reader to |4[] for recent reviews of the subject, we will start by just recalling 
the essential points needed to put the present work in the correct perspective. 

The basic postulate of PBB cosmology is, at first sight, a shocking one: our Universe 
would have originated from an "anti-big bang" state, which was essentially empty, cold, 
flat, and decoupled. The claim is that, from such innocent-looking initial conditions, a rich 
Universe can originate thanks to two distinct mechanisms: 

i) A classical gravitational instability amplifies tiny initial perturbations, inevitably push- 
ing the Universe towards a singularity in the future (to be later identified with the 
standard big bang) through a phase of accelerated expansion (inflation) and acceler- 
ated growth of the coupling. We refer to this phase as dilaton-driven inflation (DDI). 

ii) Quantum fluctuations are amplified during DDI according to the phenomenon by which 
perturbations freeze out as their wavelength is stretched beyond the Hubble radius. 
This is how we are able to produce a hot big bang at the end of DDI. 

Obviously, in order for the whole scenario to be viable, a mechanism has to be conceived 
to produce the exit from the dilaton-driven phase to the radiation-dominated phase of stan- 
dard cosmology. This is the so-called exit problem ||, on which we will have nothing new 
to say in this paper. Rather, our attention will be focused on the pre-big bang classical 
epoch, characterized by small couplings and small curvatures (in string units). The big 
simplification that occurs in this epoch is that the field equations are basically known since 
they follow from the low-energy tree-level effective action of string theory. 

One would like to discuss the most general solution to the field equations and check 
under which conditions PBB inflation takes place and is sufficiently efficient to produce 
something like the patch of the Universe presently observable to us. There have been claims 
U that this calls for highly fine-tuned initial conditions. On the other hand, arguments 
given in 0, have suggested the following interesting possibility/conjecture: the generic 
initial Universe that is able to give dilaton-driven inflation is one which, in the asymptotic 
past, converges to the Milne metric with a constant dilaton. Such an initial state does not 
look generic at first: it is so, however, in a technical sense, i.e. in that the general solution 
that develops PBB inflation is claimed to depend on as many arbitrary functions of space as 
the most general solution does. If confirmed, this would mean that PBB inflation covers a 
non- vanishing fraction of the total phase space (the space of all classical solutions) of string 
theory. 

The main purpose of this paper is to provide a non-trivial check of the above conjecture. 
We shall not address the full inhomogeneous problem, for the moment, which would require 
a much stronger numerical effort; we will limit our attention, instead, to the case of a 
spherically symmetric Universe. We thus consider the problem of how a small spherically 
symmetric lump (alternatively a shell) of energy affects the otherwise trivial evolution of 
Milne's metric. We mention that other interesting questions in spherically symmetric pre-big 
bang cosmology have been recently addressed by Barrow and Kunze |§. 
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In Section 2 we define our choice of gauge for the metric and write down the field 
equations. As usual these break up in two sets: constraints, containing only first time- 
derivatives, and (constraint-preserving) evolution equations. Luckily, we are able to solve the 
constraints in closed form and to reduce the equations to four first-order partial differential 
equations (PDEs) in one time and one space. In Section 3 we apply the general gradient 
expansion method jL0[ to our particular case and construct analytic asymptotic solutions 
near the singularity. In Section 4 we outline our technique for numerically solving the 
equations and discuss some subtleties needed to avoid possible singularities at r = 0. The 
numerical results, both for a spherical lump and for a shell, are reported in Section 5, where 
comparison and matching with the analytic asymptotic formulae derived in Section 3 are 
also made. Finally, in Section 6, we interpret our results in the string frame and discuss the 
relevance of our work to the fine-tuning question raised in ||. 



II. FIELD EQUATIONS AND ELIMINATION OF THE CONSTRAINTS 



In this paper we shall limit our attention to the graviton-dilaton system in four space- 
time dimensions. Other fields, such as the axion and various moduli, are taken to be frozen. 
The low-energy tree-level effective action in the physical (string) frame (whose metric is 
denoted by G) reads pTT] , JT2] ]: 



h- 1 S cS = ^ J d*x ^G e-* (R(G) + G» v d u 



(2.1) 



where A s , the string- length parameter, will play no role in the classical regime discussed 
here. We shall come back to its role in the later evolution in Section 6. The Einstein frame, 
which turns out to be somewhat more convenient for solving the equations, is related to the 
string frame by RitG^ = e^'^g^, with O denoting the present value of the dilaton. In 
that frame, defining 8nlp = exp(0 o )A 2 , the action becomes (12 



tr 1 S, 



off 



16tt/ 2 p 



d'x^g (R(g) - i<T d. 



>d v (/>) 



and leads to Einstein's equation (after eliminating its trace): 



and to the dilaton evolution equation: 







(2.2) 



(2.3) 



(2.4) 



For spherically symmetric cosmological solutions, depends just on time t and on a radial 
coordinate £. In the synchronous gauge we can use spatial coordinates such that the metric 
takes the form (see e.g. jl3| ): 



ds i = _ d f + e M(,t) di 2 + e 2 ^ duo 2 . 



(2.5) 



where du 2 = d0 2 + sin 2 dip 2 . Let us denote by a dot (a prime) differentiation with respect 
to t (£), introduce 
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g/g = X = 2a + A$ , A = 2 (d - /?) , (2.6) 

and recall the expression for the three-curvature on a constant t hypersurface 

®R = 2e" 2/3 - 2e~ 2a (2/?" + 3/3' 2 - 2a' (?) . (2.7) 

In terms of these, the full set of field equations consists of 
i) the Hamiltonian and momentum constraints: 



i( X 2 -A 2 ) + 2( 3 ) J R = 2 + e- 



2a j/2 



2/3A-f( x -A)' = 0'0, (2.8) 

ii) the evolution equations for the metric 

X=4(x 2 + 2A 2 )-^ 2 
A + | X A = e" 2a 0' 2 + 2e" 2/3 + 2e" 2Q (J3" - a'/3') , (2.9) 

and, finally, 

Hi) the equation for the dilaton 

+ \x4> = e~ 2a (0" + (2/3' - a'W) . (2.10) 

This last equation can be shown to be a consequence of the previous ones and we shall 
therefore ignore it. One can also show that the constraints are conserved in time, thanks to 
the explicit form of the evolution equations. The Hamiltonian constraint can be written in 
terms of the well-known Q parameter (the fraction of critical energy density) as: 

2 + e" 2a 0' 2 ^ , , 

n = P / P<r = 3 y x3 _ Aa y -1 + 6^. (2.11) 

A nice simplification, which occurs in the spherically symmetric case, is that the constraints 
Q2.8p can be solved explicitly in terms of and 0'. Given the quadratic character of the 
constraints, this leads to a four-fold ambiguity, reminiscent of the two-fold ambiguity of the 
homogeneous case, which we resolve by imposing (in accordance with the PBB postulate) a 
monotonic behaviour for at least from some time on. We then get: 

= VA 2 + B + vM 2 - B >0 

0' = e a (vM 2 + B - VA 2 - B) , (2.12) 

where the auxiliary quantities A and B are given by 

A 2 = Ux 2 ~ A 2 ) + e~ 2/3 - e~ 2a (2/3" + 3(3' 2 - 2d 0) 
B = e- a (P'A-l( X -A)'). (2.13) 

In conclusion, the final system of evolution equations contains only the metric and can 
be written in first-order form as: 



4 



« = |(X + 2A) 

/j = |(x-A) 

X = -|x 2 - |A 2 - 2(A 2 + 

A = -1 X A + 2 e" 2/3 + 2 e- 2a (/?" - c//3') + 2(A 2 - vM 4 - 5 2 ) , (2.14) 



with A, i? given in Eq. ( 2.13|) . By using the above relations and the expressions for and 



(f)' given in Eq. ( |2.12| ), it is explicitly checked that the "integrability condition" (</>)' = (</>')• 
holds. 

An essentially trivial solution to all equations is given by the so-called Milne metric, 
known in the literature [[L4| because it represents the very late time behaviour (with Q — > 0) 
of all subcritical Universes (with Q < 1). In our gauge, Milne's metric corresponds to: 

a = ln(±t) , (3 = a + lnsinh£ , cb = constant . (2-15) 

It is well known that Milne's metric can be brought to Minkowski's form through a coordinate 
transformation that maps the whole of Milne's space-time into the (interior of the) forward 
or backward light cone (depending on the sign of t) of Minkowski space-time. The conjecture 
advanced in M is that pre-big bang inflation finds its generic origin in open cosmologies that 
approach Milne as t — > — oo. This is why we will choose initial data (at some finite negative 
t) very close to Milne's Universe, by inserting some lump of energy/momentum through a 
non-trivial dilaton. 



Generic initial data in the spherically symmetric case are known [15| to depend on two 
arbitrary functions of the radial coordinate £. It is instructive to see how these two arbitrary 
functions appear in our initial data. Apparently, in order to select a definite solution to Eqs. 
( j2.14j ), one has to provide four functions of £ on a Cauchy hypersurface, i.e. a, /3,x, and 
A at some initial time. However, two of these can be eliminated by the two residual gauge 



transformations that keep us inside the synchronous gauge. We may thus take physically 
distinct initial data to be given by A and B of Eq. ( |2.13|) . Appropriate combinations of A 
and B represent the initial energy and pressure density of (dilatonic) matter. Equivalently, 
we may observe that the general spherically symmetric solution of the dilaton evolution 
equation Eq. ( |2.10|) is given (in the case of Milne's metric) by ||i~6f 



sinh£ 



y°° dp b p e ipr r^J^A + constant , r = -ln(-t/T ) (2.16) 



and thus depends on two real functions of p (since b p is complex). 

This latter observation allows us to make a more general remark: string vacua are usually 
identified with two-dimensional conformal field theories (CFTs) resulting from a set of (3- 



function equations ||11|| , fl2|l . However, within a particular CFT, we can construct marginal 
(i.e. (1, 1)) vertex operators representing physical excitations (particles) propagating in the 
given background. A physical dilaton vertex operator depends upon two (real) functions of 
space, while a physical graviton vertex depends on four. It is tempting to identify these two 
(resp. six) degrees of freedom with the classical moduli of our spherically symmetric (resp. 
generic) solutions, since we can indeed expect a correspondence between marginal operators 
and the structure of moduli space in the neighborhood of a CFT. 
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In giving initial data, two precautions have to be taken: i) we should verify that, at least 
initially, A 2 > \B\, so that all the square roots appearing in Eqs. ( 2.14|) are real. Consistency 



requires that this constraint be maintained through the evolution; ii) there is an apparent 
singularity at £ = in Puune and in its derivatives. Of course the singularity is perfectly 
canceled for Milne's metric in all equations, but one may encounter numerical problems -or 
even genuine physical ones- if one is not careful about the way to perturb Milne. We believe 
that the correct way to avoid curvature singularities and to achieve a smooth algorithm is 
to insist on an ansatz of the form: 

a = a(£ 2 ,t) , /3 = a + lnsinh£ + £ 2 5(£ 2 ,t) , (2.17) 

giving 

X = 6d + 4£ 2 5, A = r 2 A = -25. (2.18) 

It can be checked that, rewriting all equations in terms of a, 5, x, A, and their derivatives 
with respect to £ 2 , all singularities at £ = are automatically removed and that the structure 
of the ansatz is maintained during the evolution. A possible way to construct sets of "good" 
initial data along these lines is discussed in the appendix. 



III. ANALYTIC ASYMPTOTIC SOLUTIONS NEAR THE SINGULARITY 



In this section we use, as in [0], ||, the gradient expansion technique [[UJ to construct 
analytic asymptotic solutions to our spherically symmetric field equations. This exercise is 
important in several respects. Firstly, it illustrates the gradient expansion technique in a 
situation where the momentum constraints can be explicitly solved. Secondly, it will allow 
a very non-trivial check of the numerical method. Finally, it allows the numerical solutions 
to be continued in a region where the computation undergoes a critical slowdown because of 
the singularity. Hence, in the problem at hand, analytic and numerical methods very nicely 
complement each other. 

Assuming that spatial gradients become subleading near the singularity, we can simplify 
Eqs. (|2~8j), ( |2.9D to the following form: 

2 = |( X 2 -A 2 ) , 0'=(2/^-§( X -A)')/0, 

X + |x 2 = , A + | X A = 0. (3.1) 

Following the procedure of Refs. [0, we first solve the evolution equations and obtain: 

The Hamiltonian constraint can be easily solved. Choosing the branch corresponding to a 
growing dilaton we obtain: 

<K£, *) = MO - ^ v/I3 ^ ln (*/*o(0 - 1) , t < to • (3-3) 
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We now integrate once more the evolution equations to obtain the metric: 



« = 7(0 + |(l + 2A(0)ln(t/t (0-l) , 

= *(O + !(l-A(fl)ln(f/to(0-l) • ( 3 ' 4 ) 

Finally, we solve the momentum constraint. Most terms automatically match and we 
jet the single condition: 

2v ' ] (X'/3 + XS') . (3.5) 



The general solution thus appears to depend upon four functions of space, i.e. A, 7, 5, and 
to, giving the line-element: 



ds 2 = -di 2 + (t/to - 1) 2/3(1 - A) [(t/t - 1) ZA e 2 ~< d( 2 + e 2d du 2 \ (3.6) 

Obviously, we can reabsorb 7 in a redefinition of the £ coordinate, dr = e 7 <i£. Similarly, 
the choice of equal-time slices allows us to remove the space dependence of t . We are thus 
left, as in the previous section, with just two physically meaningful functions of £ for the 
characterization of our dynamical system, as should be the case for the general spherically 
symmetric solution [15[ . 



2A ^7 J f 2 , J26, 



Actually, the situation is a bit more subtle as far as the £ dependence of to is concerned. 
Inserting the above asymptotic solution back into the exact equations, one finds that spatial 
derivatives are only subleading for sufficiently small values of t (£)- In other words, at least 
in the Einstein frame, one has to choose appropriate time slices so that to is constant and 
only then can one neglect spatial derivatives. Constructing these privileged time slices is 
not a simple problem and therefore, in this paper, we will only check (see Section 5) the 
asymptotic formulae in the vicinity of the minimum of to (i-e. near the point where the 
singularity is first reached). In Section 6 we will argue that an alternative solution to this 
problem consists in going over to the synchronous gauge in the string frame. In this case, 
spatial derivatives (at fixed string-frame time) turn out to be always subleading near the 
singularity. This suggests that it would be desirable to numerically solve the equations 
directly in the string frame. 



IV. THE NUMERICAL APPROACH 

We shall now briefly describe the essential aspects of the numerical algorithm and its 
practical implementation. The system given by Eqs. (2.14) is similar to a Hamiltonian one 



with an infinite number of degrees of freedom. Its numerical implementation must necessarily 
introduce a limit on the number of degrees of freedom, which can be done in several ways. 
Our present choice is to introduce a finite grid in the variable £ and to define the derivatives 
with respect to £ by a spectral method, based on the Fourier transform, which allows us to 
reach a high precision. Alternative techniques, such as using a symbolic language or applying 
finite elements techniques, will be considered in the near future. In our approach, a special 
care must be devoted to boundary conditions, since the Fourier transform preferably works 
in periodic boundary conditions, which are not adequate to our problem. As we explain later 
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on, this problem is solved by letting £ extend to a symmetric interval and by continuing the 
fields to the negative £-axis. Time evolution is generated via a standard integration routine 
with variable step and local error control. 

To prepare the equations to be implemented in numerical terms, it is convenient to 
introduce new fields, which describe the metric as a perturbation from a Milne background: 

a s = a + t 

P s = (3 + r — In sinh £ 

Xs = e~ T x 

A s = e- T A (4.1) 

where r is defined by t = — exp(— r). The advantage of working in the Milne background is 
that the shifted fields, unlike the background, vanish at £ = oo and we can thus introduce 
a finite volume cutoff. The boundary conditions at £ = are crucial; as we have already 
noticed, regularity at £ = is achieved if we assume Eq. (|2.17f) , in the equivalent form 



/? s (£,T) = a s (£,T)+£77(£,T) 

A,(£,r)=£S(£,r) (4.2) 

with boundary conditions at £ = 

r/(0,r) = S(0,r)=0, (4.3) 
a.(0,r) / = x.(0,r) / = 0. (4.4) 

In order to satisfy these boundary conditions and still apply the Fourier transform to com- 
pute field derivatives, we continue the fields to a symmetric interval — L < £ < +L; we 
enforce the correct boundary conditions by extending a s ,x.s to symmetric functions and 



rj, X to antisymmetric ones. In terms of the new fields, the system ( |2.14j ) becomes 



dr 
drj 

dr 
dXs 
dr 



- Y _ ly 2 _ lA 2 - A , 



— = -E - \xs^ + 2C 1 e" 2Qs f g . nh2e + - a' s (fl + coth£)J + C* A_, (4.5) 
where 



A± = 2 ( At ± a/A! - B 



6 ~ 2aS [ sinh 2 £ ~ 3 " 2/? " ~ ~ 2{3P ' S ~ a ' s) C ° the + 2a '^' s ) 
B s = e -* (A. (fl + cothO - |(y/ s - A'j) . (4.6) 
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and Eq. ( |4.2|) is understood. 

Notice that the system is identically satisfied by taking a s = 77 = E = and Xs — — 6, 
which represents Milne's background. The initial conditions we want to examine are given 
by Eq. (|A1|) (with ( = 0), which in terms of the shifted fields reads as follows: 

^0)^ln(l- cosh[ ^_ eo)] ) + (^ ^-e ) 

Xs (£,0) = -6exp{-2a s (£,0)} (4.7) 
77(£,0) = £(£,0) = 0, 

where /i is a suitable constant. We have implemented this system of partial differential 
equations in MATLAB which offers very efficient built-in routines for integrating ordi- 
nary differential equations, a fast Fourier transform and an integrated graphic environment .[] 
The only tricky point in the numerical treatment of these equations refers to the delicate 
cancelation mechanism that makes the solutions regular at £ = in spite of the apparent 
singular nature of the terms involving coth£ and sinh _2 £. This cancelation, which takes 
place in exact arithmetic, is spoiled by finite precision arithmetic and makes the solution 
singular after a short time evolution. The remedy that we adopt consists in enforcing regu- 
larity near the origin by a polynomial interpolation and a Fourier filtering, which truncates 
high frequency components above a certain cutoff. A consistency check is given by solving 
the equations with various initial ansatze corresponding to £0 — or £ > 0. The former case 
represents a spherically symmetric perturbation concentrated at the origin, while the latter 
spreads out over a spherical shell and naturally avoids the singularity. A similar behaviour 
is indeed observed in both cases with fields a and x blowing up as Q approaches 1. 

V. NUMERICAL RESULTS 

We report some preliminary results obtained by running our code with the initial ansatz 
given by Eq. ( |4.7| ) and for several values of the parameters (£o>£> /-0; a typical run involves 
a grid with N = 4096 points and takes less than an hour on a modern workstation. We set 
a finite-size L « 10 /e + £0 and an ultraviolet cutoff at half way to the size of the Brillouin 
zone. The routine ode45 of MATLAB can be used with its standard setup; a special care must 
be devoted to the choice of the absolute tolerance parameter which should be set differently 
for the various fields which have very different scales. The evolution starts at negative t 
from the ansatz of Eq. ([4.7]) ; we report both cases (£0 = and £0 > 0), which indeed show 
a similar pattern (see Figs.(jl],0)). In Figs.(0,fD both fl and the dilaton are shown. 



1 A copy of the code is available at the URL tittp : //www. f is .unipr . it/~onof ri 
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We examined the numerical solution in order to check the validity of the approximation 



introduced in Sec.|T|. In particular we were able to check the asymptotic behaviour in the 
region where Q is very close to 1, where Eqs. ( |3.2| ) and ( p.4[ ) predict a linear regime in time 
for the following expressions 



W 1 = l/x 

W 2 = exp{a + 2(3} 

W 3 = exp{-(|0 + a + 2/3)/(V3 - 1)} 
(notice that the value of A is negligible). This is clearly displayed in Figj5|. 



(5.1) 



x 10" 



a) 





-0.32 



-0.3 



-0.28 



t 



FIG. 5. The asymptotic behaviour 0(io(£) — t) for the variables W n at the O peak (a) £q = 0, 
b) = 10) 



From a linear fit to the data, it is then possible to extract information about the unknown 
functions to(0 an d A(£). We find that, as expected, to(£) has an extremum near the fl 
peak with a small curvature £(£)", which is responsible for the small deviations from the 
asymptotic estimates of Eq. (|3.2|) (see Fig.^). 

In all the cases we have examined so far, the function A(£) turns out to be very small (see 
Fig.|7|). but this is probably due to our initial starting point with A = 0. We are currently 
exploring a wider class of initial data, in order to find the generic properties of the solutions. 
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FIG. 7. The function A(£), (oj Co = 0, b) £ = 10). 



Finally, our code can also be run backwards i.e. towards t = — oo. In this case, however, 
we encounter a problem: although keeps decreasing in magnitude as t becomes more and 
more negative, it also starts to oscillate, as expected ||. As a result, the constraints cannot 
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be always solved with the sign determination given in ( |2.12|) . We have been able to solve this 
problem by imposing the constraints only on the initial data and by otherwise working in 
an enlarged phase space containing also <fi and <fi. This procedure has the further advantage 
of allowing a non-trivial check of the numerical precision by verifying that the constraints 
are conserved. 

We expect to present results obtained by following this better procedure in the near 
future. 



VI. DISCUSSION 

In order to discuss better the physical meaning of our result it would be convenient to 
transform them back to the original string frame. It is difficult to make the change of frame 
numerically and, indeed, we think it would be easier to start directly the whole exercise in 
the new frame. Fortunately, we can achieve a semi-quantitative understanding of the string 
frame solutions by noticing that: 

i) initially the two frames coincide, since the dilaton is constant; 

ii) our numerical results, as the singularity is approached, fit very well with the analytic 
asymptotic formulae given in Section 3. 

The latter observation allows us to perform the passage to the string frame analytically (and 
of course approximately) near the singularity. Following Ref. we find rather easily that 
the asymptotic metric and dilaton in the string-frame synchronous gauge are given by 

ds 2 = -dt 2 + [(t/to - l) 2a « e 2 ^ d£ 2 + (t/to - lf ae ' v e 25 du 2 ] , 
<K£, t) = <M0 + 7<f> In (t/to(0 - 1) , * h (6.1) 

where, for notational simplicity, we kept denoting by t the string frame time and by io(0 the 
location of the singularity. The important exponents a% and 7^ are completely determined 
by the function A(£) appearing in Eq. Q3.2p . We find: 



1 2x/r^P- y^A 2 -2A (V3 + y^A^) 
73 2 + A 2 

1 2 V / T^A 2 - v^ + A {V3 + VT^W) 



! ~ v/3 2 + A 2 

2(1 - A 2 ) + 2 v / 3yr^A 2 
2 + A 2 



^ = ~~ o ■ x2 • (6-2) 



It is reassuring to check that the oti satisfy Kasner's condition: 

£of = af + 2a^ = l (6.3) 

and that, furthermore, 

7 *-$>i = -l (6.4) 
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as in the homogeneous case. Also note that, in the limit A < 1, aj -> — l/v3. This 
corresponds to the isotropic limit. From the above formulae we see that, in the string frame, 
the metric is typically (super) inflationary since the inequalities — 1 < on < are almost 
always fulfilled. It is easy to check that, as a consequence, the string frame asymptotic 
solution is reliable even if one does not choose time slices corresponding to to ~ constant. 
The expansion rate is maximal for the regions where A is very close to zero, i.e. for very 
isotropic regions. 

We conclude that, as anticipated in refs. @, [§], regions with the correct dilaton per- 
turbation do undergo a superinflationary expansion, become asymptotically flat and, most 
likely, isotropic. It would of course be most interesting to attempt to generalize our com- 
putations to the case of several lumps and/or shells located randomly in space and with 
random initial parameters and see what kind of chaotic inflationary Universe will result. 
This, unfortunately, appears to demand going beyond the spherically symmetric situation 
we have considered in this paper. 

In conclusion, our numerical approach appears to confirm that, at least in the case of 
spherical symmetry, dilaton-driven inflation naturally emerges as a classical gravitational 
instability of small perturbations around the trivial (Milne) vacuum. Nonetheless, before 
being able to address/answer completely the criticism expressed in ||, the present work 
needs to be expanded in at least two directions: i) classical phase space has to be swept more 
systematically, in particular away from the case of spherical symmetry; ii) the behaviour at 
very early times has to be thoroughly investigated so that the possibility (see N. Kaloper et 
al. 0]) that quantum fluctuations at very early times may destroy the homogeneity needed 
for turning on dilaton-driven inflation is properly assessed. As discussed at the end of 
Section V, this latter problem requires a new way of enforcing the constraints on which we 
are presently working and hope to report soon. 
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APPENDIX A: A PROCEDURE TO GENERATE A CLASS OF ADMISSIBLE 

INITIAL DATA 

In this appendix we illustrate our procedure for choosing the small initial perturbation 
of Milne's metric in such a way as to fulfill the necessary positivity constraints as well as 
regularity at £ = 0. The physical idea is quite simple: since we know that Milne's solution is 
regular and satisfies all the constraints, we construct initial data from a small deformation 
of Milne. Notice that there is no need that the deformation satisfies the evolution equations. 

Consider the following deformation of Milne's metric: 
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(3 = a + In sinh ^ + 



T 2t2 



t 2 coshe(£ — £1 



(Al) 



Such a deformation is localized at £ = £o, with some anisotropy concentrated at £1. The 
deformation rapidly drops to zero far from these two values of £. The symmetrization is 
needed in order to bring the ansatz in the form discussed ion Section 3. We now define the 
initial data from Eqs. ( |Al| ) and their first time derivatives at the initial time ti <T . 

We have performed the computation of the initial data using Mathematica and checked 
that the positivity and regularity constraints are identically satisfied in £ for a wide choice 
of the parameters U/Tq, e, £ , £i, (, e.g. for the sets used in section 5: 

U/T = 10, e = 0.2, £ = 0, C = 0; 
U/T = 10, e = 0.2, £ = 10, C = . (A2) 

Indeed, positivity is most sensitive to the parameter e which cannot be much larger than 0.5. 
In the limit U/Tq — ► oo the formulae for the initial data simplify considerably and positivity 
constraints can be solved in many cases analytically. 

The above initial choice, that we have considered for our numerical study, is a special 
case of a class of admissible initial data (we call a given choice of initial data admissible if 
it satisfies the constraint A 2 > \B\) which can be characterised as follows: 
theorem: Let /(£) be twice continuously diSerentiable with 

(i) /(0 > fo > i; 

(ii) f (0) = 0; 

(hi) f"(0 < e 2 /(0 for some e < 1. 
Then the initial data given by 

a s = | ln(l - /(C)- 1 ) , X s = -6 e~ 2a ° , V = E = . (A3) 

are admissible if 

/-i < max[ 5 ~ 4£ ~ 2£ ' , 1 - ( — -) 2 ] . (A4) 

Jo ~ L 6-4 £ -2e 2 ' K 3-2e-e 2 ' 1 V ' 

We shall first prove the following 

LEMMA: Let /(£) satisfy the conditions of the previous theorem. Then 

f'(£) 

-£<^<£tanh £ £ (A5) 



proof: Let = /'//; then it holds 



<e 2 -g{if. (A6) 



First of all let us show that if g(£) < —e for some £, g would develop a singularity at a finite 
£: under the assumption and the previous inequality Eq. (|A6|) it would follow 



16 



g'(0 



> 1 



which can be integrated to give 



9(0 dg 



9(0 9 2 ~ £ 2 



(A7) 



(A8) 



Since the integral is convergent, g — > — oo at a finite £; / being regular and > 1 this is a 
contradiction. Assume next that > etanhe^ at some £; by continuity we may assume 
< e. It follows 

2 g ^L < i ; (A9) 
e -0(f) 

by integrating from £ = we find tanh -1 (g(£)/£) < e£, a contradiction. This completes the 
proof of the lemma. 

To prove the main theorem, let us insert Eq. ( |A3| ) into the expression of A 2 S ±B S (Eq. (|4.6|) ); 

we get 



(/ - 1) exp{2a s } (A\ ± B s ) 



(A10) 



f" 2/ -5/4 f'Y f ,,,,„/'// 

3 -J + ^r{7) - 2 7 cothe±2 7V7^i- 

We can apply the lemma at once and see that the l.h.s. is bound from below as follows 

2 



l.h.s. >3-e-2e+2 



r 



± 2 



r f 



jj /v/-i 

or, we may complete the square before applying the inequality, 



l.h.s. > 3-e z 



*°- k A + \7^A. 



(All) 



(A12) 



and the result follows. Notice that the theorem covers the cases that we have considered in 



section IV; the domain of parameters for which the data are admissible is shown in Fig| 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 8. Admissible values for fi, e. 
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The class of initial data which has been defined here depends on a single arbitrary 
function. The general expression of / can be easily derived by setting e 2 f — /" = p(£) > 0. 
The solution is clearly 

/(0 = /i- 1 ooBh(eO - f e~ l sinh (e(£ - £')) p(0 d? . (A13) 

Jo 

If / °°p(O rf^ is small enough, then we get admissible initial data. As stressed in Sec.||, it 
would be interesting to have a general characterisation of suitable Cauchy data. 
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